Structure of fermion nodes and nodal cells 



o 
o 

(N 

c 

a 

o 

(N 



Lubos Mitas 1 

'Department of Physics and CHiPS, North Carolina State University, Raleigh, NC 27695 

(Dated: February 6, 2008) 

We study nodes of fermionic ground state wave functions. For 2D and higher we prove that 
spin-polarized, noninteracting fermions in a harmonic well have two nodal cells for arbitrary system 
size. The result extends to other noninteracting/mean-field models such as fermions on a sphere, 
in a periodic box or in Hartree-Fock atomic states. Spin-unpolarized noninteracting states have 
multiple nodal cells, however, interactions and many-body correlations generally relax the multiple 
cells to the minimal number of two. With some conditions, this is proved for interacting 2D and 
higher dimensions harmonic fermion systems of arbitrary size using the Bardeen-Cooper-Schrieffer 
variational wave function. Implications and extent of these results are briefly discussed. 

PACS numbers: 02.70.Ss, 03.65.Ge 
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A fermion node is a subspace of fermion configurations 
for which a real wave function describing the fermionic 
system vanishes due to the antisymmetry. In general, for 
N spin-polarized fermions in d dimensions the fermion 
node is a (dN — l)-dimensional hypersurface given by an 
implicit equation ^f(R) — 0, where $ is the wave function 
with fermion coordinates R = (ri, tn). The location 
of the nodal manifold is of key importance for quantum 
Monte Carlo (QMC) methods [HUSH since the exact 
node enables us to solve the stationary Schrodinger equa- 
tion with computer time scaling as a low-order polyno- 
mial in N. Remarkably, even rather crude Hartree-Fock 
(HF) or post-HF wave function nodes, routinely used in 
the fixed-node approximation QMC methods, provide 
95% of the correlation energy in a variety of systems 
with hundreds of valence electrons 0,0, 0-Q However, 
to reach beyond this level of accuracy has proved to be 
challenging because of limited understanding of fermion 
nodes 0, 0, B B ■ 

The exact nodes for interacting systems are known 
only for a very few two-electron triplet atomic states 
|.4ll(l| and, very recently, the exact node of three-electron 
4 S(p 3 ) state has been discovered ^lj. ^ turned out that 
the nodes of these high-symmetry states have rather sim- 
ple topologies and divide the configuration space into two 
compact nodal cells in which the wave function is posi- 
tive or negative ("plus" or "minus" cell). The two nodal 
cells were found also in noninteracting spin-polarized 2D 
and 3-D fermions with up to 200 particles using a nu- 
merical proof The same noninteracting but spin- 
unpolarized systems trivially have four nodal cells since 
the wave function is a product of spin-up and spin-down 
determinants. Interestingly, analysis of interacting unpo- 
larized few-particle systems 
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has revealed that 

the electron correlation can change the node topologies 
and number of nodal cells. For example, the correlated 
wave functions of the Be atom or N2 molecule exhibit two 
nodal cells while the corresponding Hartree-Fock ones 
have four. It was therefore conjectured that the bisection 
of the configuration space into the two nodal cells might 



be a generic property of fermionic ground states (with ex- 
ceptions discussed later) . In this Letter we prove that for 
2D spin-polarized noninteracting harmonic fermions the 
ground state node divides the configuration space into 
the minimal number of two nodal cells for any system 
size. The proof method extends to higher dimensions and 
carries over to other models such as fermions on a sphere, 
in a periodic box, for atomic states, etc. We show that 
the same holds, in general, for interacting systems using 
correlated Bardeen-Cooper-Schrieffer wave functions and 
we briefly discuss the implications. 

We recall two basic properties of fermion nodes derived 
by Ceperley 0. a) Nonde generate ground state wave 
functions satisfy the so-called tiling property which states 
that by applying all possible particle permutations to an 
arbitrary nodal cell one covers the entire configuration 
space. Note that this does not specify the number of 
nodal cells, b) Consider three particles k in a spin- 
polarized system with wave function ^(i?). We call the 
particles i, j, k connected, if there exists a triple exchange 
ijk — ► jki path that does not cross the node, ie, > 
along the exchange path. More connected particles, 
such as the following six ones, can form a cluster: <J*J . If 
there exists a point R t such that triple exchanges connect 
all the particles into a single cluster then ^S(R) has only 
two nodal cells. The tiling property implies that once 
the particles are connected for R t the same applies to 
any point in the cell as further explained in Ref. 

We illustrate the properties of nodes in ID using 
spin-polarized fermions in a ID harmonic oscillator well. 
The wave function is a Slater determinant N) — 

det{0 fc (r 2 )} = -An i e-^ /2 det{l,2a:,... > fr w _ 1 (a;)} where 
H n (x) is a Hermite polynomial of degree n and A is the 
normalization. We omit the prefactors and transform the 
Slater matrix to monomials so that the wave function is 
given by the Vandermonde determinant 

* 1£) (1, ...,JV) = det{l,x,x 2 ,...,x N - 1 } =]J{x j -x % ) 



(1) 
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The node is encountered whenever two fermions pass 
through each other and the wave function has N\ nodal 
cells since any permutation requires at least one node- 
crossing. In general, the derived node is exact for other 
ID models including systems with interactions. 

Now consider spin-polarized fermions in a 2D har- 
monic well. The one-particle states are simply (f> nm = 
C nm H n (x)H m (y) , n, m = 0, 1, ... where C nm includes the 
gaussian and normalization which are absorbed into a 
common prefactor and omitted. The Slater matrix ele- 
ments can be rearranged to monomials and we write 



For example, for M = 2 the wave function reads 



V 2D (l,...,N) = det{l,ar,y, 



,x n y m , 



•} 



(2) 



The closed-shell states and the system size are labeled 
by M — 1, 2, ... where n + m < M with the number of 
fermions given by N = (M + 1)(M + 2)/2. 

Using induction we prove that the wave function in 
Eq. has only two nodal cells for any M > 0. This 
is indeed true for M — 1 with three-particle wave func- 
tion 4 , 2_d(1,2,3) = det{l,x,y}. In order to show this, 
it is convenient to extend the particle 2D coordinates 
by a "dummy" third dimension as = (xi,?/j,0). Then 
^20(1,2,3) = z • (r 2 i x r 3 i) where z is the unit vec- 
tor in the third dimension and r.^ = — rj. Clearly, 
there are only two nodal cells since the set of vectors 
Zo, r2i, is either left- or right-handed, and the node is 
encountered whenever the three particles are collinear, ie, 
F2i x r3i = 0. Also, the particles are connected by triple 
exchanges without node-crossing (eg, rotate an equilat- 
eral triangle). 
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where the particle positions are given in Fig^ (our point 
Rt). Clearly, all the last-row elements containing yo can 
be eliminated by adding a multiple of the third row. In 
a similar way we eliminate all the matrix elements con- 
taining yo, and the determinant factorizes as 

2 

*2D(i,...,6)=n^-^) n! * ic ( i ' 2 ' 4 )* 2o ( 3 ' 5 ' 6 ) 

1=1 

(4) 

where m is the number of particles on line y — yi. Note 
that one of the factors is the 2D wave function for the 
system with M reduced by one. The same factorization 
structure is obtained for an arbitrary M and, obviously, 
it can be applied recursively. The wave function for a 
general size M is then given by 

M 

#223(1, ...,N) = ± Y[( yi - yo) ni *lD(I )*2D(l, N/I ) 
i>0 




x O x 1 x 2 ' ' ' X M X M+1 

FIG. 1: Particles arranged in a Pascal-like triangle pattern. 
The circles denote the particle positions and inside the circles 
are particle labels. 

We now consider a general system with M > 1 and ar- 
range the particles into a Pascal-like triangle pattern on 
a rectangular mesh as shown in Fig^ For this arrange- 
ment the determinant can be explicitly evaluated for any 
M by subsequent factorization of lines of particles |13| . 
Due to the space constraints and to simplify the notation 
we illustrate the factorization on a few-particle example; 
the generalization to an arbitrary size is straightforward. 



M-l 



1=0 



Hiyi-y^mili) 



i>l 



(5) 



where // 



■(I) 

l M+l-l 



denotes the labels of particles 



lying on the line I while (1, ...,N/Iq) means that the la- 
bels in I are omitted from 1, N. The sign depends on 
the number of row exchanges and on the actual ordering 
of particles. Note the translation invariance of Eq. [5j 

For the induction step we assume that triple exchanges 
connect all the particles in the system of size M. Let us 
show that the same is true for the system of size M + 1 
(see Fig^l . We first factorize out the line y = yo as given 
by Eq. [S] One of the factors is the 2D wave function for 
the system of size M, which contains all additional par- 
ticles except the one with coordinates [xm+i, yo]- By the 
original assumption, such a system has all the particles 
connected. If instead of the horizontal line we factorize 
out the vertical line (x — xq) we see that the particle at 
[xm+1i yo] is also connected, thus concluding the proof. 

The proof can be generalized in several ways. The same 
arguments apply to 3D since particles can be arranged 
into an appropriate 3D pattern and the wave function 
evaluated by recursive factorization of planes and lines. 
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In fact, the result holds and the two nodal cell property is 
correct for arbitrary dimension d > 1! The proof applies 
also to other noninteracting or mean-field models with 
polynomial entries in the Slater matrix such as fermions 
on a sphere or in a periodic box with details given else- 
where ■ The proof can be modified at least for some 
open-shell states as well while takinginto account node 
ambiguities if there are degeneracies [j, ||| . 

The proof method can be further combined with sym- 
metries for cases when a complete factorization is not 
obvious, such as for multi-shell atomic states. For ex- 
ample, spin-polarized ground states Is2s2p n ,n =1—3 
(and beyond) also have two nodal regions. We show this 
for the 6 S(ls2s3p 3 ) state in the noninteracting and HF 
limits. The wave function is written as \& a t(l, 5) = 
det{p* ls (r) 7 p* s (r),x,y,z} where p* ls (r) = pu{r) / p 2p {r) 
and p 2s {r) = P2s(r) / P2 P (r) since nonnegative p2 P (r) 
is factorized out. The dimensionless coordinates are 
rescaled by the atomic number Z and the Bohr radius 
ao as r <— Zr/ao- Let us place particle 1 at the ori- 
gin and particles 2-5 on the surface of a sphere with the 
radius 770 equal to the radial node of p2s(j~) orbital, ie, 
P2s(i]o) = 0. For such configurations we obtain 

* at (l,...,5) =/>L( J ?o)pL(°) r 32 ' (>42 x r 52 ) 

so that any three-particle exchange from 2,3,4,5 easily 
avoids node-crossing by appropriate positioning and ro- 
tations on the sphere. The particle 1 is connected by the 
exchange 123 — > 312 parametrized as i"i(t) = r]o[t,0, 0]; 
r 2 (t) = T] O [c(t),s(t),0] and r 3 (t) = r) [0,l - t,0] where 
t = (t = 1) corresponds to the beginning (end) point 
of the exchange path while c(t) = cos(7rf/2) and s(t) = 
sin(7rf/2). Setting r 4 = [0,0, 770] and r 5 = [0,0, —770] we 
find that ^S> a t is proportional to 

p* 2s (t m )c(t)(l - t) + p* s [(l - t) m ]s(t)t > 

The inequality holds for the whole path < t < 1 since 
p^si^Vo) > for < t < 1 for both noninteracting and 
HF cases. The proof can be further extended to more 
shells such as Is2s2p 3 3s3p 3 and Is2s2p 3 3s3p 3 3d 5 Q. 

Spin-unpolarized systems. Due to the product of spin- 
up and -down determinants, the number of nodal cells in 
noninteracting unpolarized systems is twice the number 
of cells of the half-filled spin-polarized counterparts. The 
proof for atomic states above then implies that the HF 
wave functions for atoms with Z > 3 up to the third-row 
elements have four nodal cells. 

Interactions in spin-polarized systems. In general, the 
shape and topology of the nodal manifold is influenced 
by interactions and many-body effects. Consider the 
lowest atomic quartet of S symmetry and even parity 
4 S e (ls2s3s) (clearly not the lowest quartet, which is 
odd 4 P°(ls2s2p)). The noninteracting wave function 
is given by * Qt (l,2,3) = det{p ls (r), p 2s (r), p 3s (r)} and 
has six nodal cells. Since 4' t(l,2,3) depends only on 



distances it is quasi-11? and the nodes are the same 
as the ones given by Eq. ^ F°r the interacting sys- 
tem the correlation is included by adding to the wave 
function the lowest (and dominant) double excitation 
2s3s — > 2p x 3p x + 2p y 3py + 2p z 3p z with a weight w. 
The correlated wave function then allows for exchanges 
without node-crossing. Define r a (t) = [0, c(i), s(t)], 
r 5 (t) = [c(t),«(t),0], r ab (t) = [1 - g(t)]r b (t) - g(t)r a (t), 
r c {t) = r ab (t)/\r ab (t)\ 7 and r d (t) = r c (i) x r (t), where 
g(t) = 3t(l - t). The exchange path 123 -> 231 is 
then r 3 (t) = r 3 (t)r a (t), r 2 (i) = r 2 {t)r c (t) and ri(t) = 
r i(t) r d(t)/\rd(t)\. The radial parts are given by 7*1 (t) = 
r] m + q(2t-l), r 2 (t) = <q m + q{l-t) and r 3 (t) = r) m -qt. 
r\ m is the mean value of the radial node of p2 S and the first 
radial node of p 3s orbitals while < q < ao\w\ (for the 
Coulomb e-e interaction w is rj —0.05). The path is or- 
chestrated so that in the region where the noninteracting 
component vanishes the correlation dominates and the 
particles become connected. This illustrates two points: 
imposing symmetries at the mean-field level can generate 
multiple cells and, in general, for d > 1 the interactions 
lift this "nodal cell degeneracy" and relax multiple cells 
to the minimal two. 

Interactions in spin-unpolarized systems. The change 
from four to two nodal cells due to interactions has 
been demonstrated for the first time on the Be atom 
using the two-configuration correlated wave function 
and the connected cluster construction adapted to spin- 
unpolarized systems. Consider a simultaneous exchange 
of one (or more) pair(s) of spin-up particles and one (or 
more) pair(s) of spin-down particles. For noninteracting 
wave functions such simultaneous pair exchanges imply 
that the node will be crossed once or multiple times. If 
there exists a point Rf such that during the simultaneous 
spin-up and -down pair exchanges the inequality | ^ | > 
holds along the whole path, then the wave function has 
only two nodal cells. 

Consider a singlet state of 2N particles in a 2D har- 
monic well with particles interacting by pair potential. 
With some restrictions, we show that the correlation 
included in the Bardeen-Cooper-Schrieffer (BCS) pair- 
ing wave function 0, E| given by ^scs(l) — , 2JV) = 
det[<f>(i, j)] is enough to eliminate the noninteracting 
four nodal cells and fuse them into the minimal two. 
Here $(i,j) = i) is a singlet pair orbital for if 
and j I fermions and we decompose it into noninter- 
acting and correlated components <fr(i,j) = $>o(i,j) + 
&corr{i, j)- Using one-particle orbitals we can write 

j) = S^ TO =0^nm( i )^nm(j) while $ corr (i,j) = 

E„ +m >A/ c nm (f>nm{'i)(f>nm(j) where {c nm } are variational 
parameters. 

We illustrate this first on a six-particle singlet, M — 1, 
2N = 6. In cylindric coordinates (r, tp) we find $o(i,j) — 
1 + 2r^j cos ifij where <fij = <pi — <fij, with the gaussians 
and normalization omitted. $ C orr(*,j) is constructed as 
a sum of orbitals from the next unoccupied shell (n+m = 
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FIG. 2: (a) Positions of six 2D harmonic fermions with two 
particles at the origin and two pairs on circles with radii r a ,rb- 
(b) Positions of the spin-up particles for the M = 2, 2N = 12 
particle singlet. The spin-down positions pattern is similar 
except for the different radii and angles. 



2) and we find $ cor r(i, j) = a[2{rir cosipij) 2 - rf - rf\, 
where a variational parameter a has been included. We 
place the particles as in Fig^a). For such a configuration 



BCS 



1, 



,6) = 8ar a r b cos</5[2(r a r fc cos</3) 2 



The rotation of the system by 7r exchanges the pair of 
particles for each spin, and since the wave function is 
rotationally invariant, it is enough to show that it is 
nonzero for a single point, eg, r a = Tb = 1, <~p- = 7r/4. 
For a = the wave function vanishes since the particles 
lie on the noninteracting node. However, for any inter- 
action which does not break the rotation invariance and 
gives the BCS wave function has only two nodal 

cells. 

Remarkably, this can be generalized to an arbitrary 
size. Assume a closed-shell singlet with the total num- 
ber of particles 2N = (M + 1)(M + 2) where N is even 
(for N odd the derivations are the same after placing 
one particle of each spin to the origin, see Fig|2fa)). 
We form N/2 pairs of particles in each spin subspace 
so that, say, the pair i f, (i + 1) t has coordinates given 
^y(n><Pi) = (rk,<Pk), (n+i,<Pi+i) = (rk,<Pk+n), see Fig. 
Hfb). Here k = 1, N/2 labels the pairs in the spin-up 
channel; the spin-down particles are placed similarly and 
labeled by the pair index I = N/2 + 1, N. In this con- 
figuration the particles lie on the noninteracting node 

since det{$ («,j)} = <iet{J2 n+m <M <t>nm(i)<l>nm(j)} = 
det{4> nm (i)}det{ip nrn (j)} and the rotation of the sys- 
tem by 7r crosses the nodes in both spin channels so 
that both Slater determinants vanish due to the rota- 
tion invariance. Now, if all the N pair distances and an- 
gles rk,ri, ipk, tpi are distinct, then each of the matrices 
{4>nm('i)}, {4>nm(j)} has exactly one linearly dependent 
row, ie, their ranks are N — 1. This can be verified di- 
rectly for small values of M and then using induction for 
any M. Consequently, the matrix {<£>o(*, j)} has linear 
dependence in one row and one column, ie, it has the 
rank of N — 1 as well. In general, adding virtual states 
through ^corr(hj) provides independent rows/columns 



(eg, M + 1 independent rows/columns from the first un- 
occupied shell) which eliminate linear dependency so that 
det{$ (i, j) + &corr{i,j)} is nonzero. Assuming that the 
interactions do not break the rotation invariance, the cor- 
related BCS wave functions have only two nodal cells. In 
fact, this can hold even if the invariance is broken but 
one would need to show it for the entire exchange path, 
not only for a single point. The proof extends to d > 2 
and to other models as well, with details given elsewhere 

Q 

For the classes of fermion systems studied in this work 
the two nodal cell property indeed appears as a generic 
feature. That brings up an interesting question: when 
might this property not apply? We mention just some of 
the possibilities: i) additional symmetries and/or bound- 
ary conditions can generate additional nodal cells; ii) 
nonlocal or very strong/singular interactions can reorder 
the states (eg, an excited state becoming the ground 
state) or significantly change the nodes; iii) for systems 
with strong correlations, large degeneracies at the Fermi 
level, or at quantum phase transitions the properties of 
nodes are far from clear and require further study. 

We believe the presented analysis and proofs provide a 
significant step forward in our understanding of topolog- 
cal properties of fermionic wave functions. The results 
offer an elegant and unifying framework for several pre- 
viously conjectured or numerically investigated features 
and apply to arbitrary size for both noninteracting and 
interacting systems. The presented techniques extend to 
other models and open exciting perspectives for studies 
of many-body effects which are currently out of reach. 
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